A single measurement of a quantum many-body system of bosons 
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Here I propose an approximate way of simulating the outcomes of a single- experiment density 
measurement that is performed on a state of JV bosons. The approximation is accurate if occupation 
of single-particle modes is macroscopic. 
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INTRODUCTION 

In order to simulate a single measurement of a position 
of a quantum particle described by a wave function <j>(x) 
it is enough to randomly draw a position x = £ with the 
probability density |0(x)| 2 . The outcome of the measure- 
ment, the particle at a point x — £, is different from the 
statistical average over many such measurements, given 
by |<^(x)| 2 . This kind of difference, trivial in the case 
of measurements performed on a single particle, becomes 
especially interesting in the case of quantum systems that 
consist of many bosons. One prominent example is given 
by two colliding Bose-Einstein condensates 0] , where ev- 
ery single measurement of the system reveals an interfer- 
ence pattern while no pattern is present in the average 
over many measurements. 

For a generic quantum state of a many-body system it 
is inexpensive to come up with a prediction for an aver- 
age outcome of many measurements. On the other hand, 
finding out possible results of a single measurement is 
extremely difficult. If the state is given by a many-body 
wave function 4>{x\, ...,xn), where Xi stands for a coordi- 
nate of an i-th particle, one needs to draw a set of N po- 
sitions (one for each particle) from the iV-dimensional[3 
probability density \(fi(xi, ...,Xn)\ . For large number of 
particles, A, the direct sampling of the corresponding 
multidimensional probability density is very difficult, if 
possible at all. The direct sampling of the iV-dimensional 
probability density can be replaced by sampling of N one- 
dimensional conditional probability densities |2| ■ In prac- 
tical applications, however, this clever method is suitable 
to handle large number of particles only if they occupy 
very few modes. 

In the following I propose an approximate method of 
simulating outcomes of single measurements. It is de- 
signed for those quantum many-body states that involve 
many macroscopically occupied modes. 

A reader, novice to the subject, may establish a neces- 
sary background by contemplating two first pages of Q 
and three first pages of [2| ■ 

Suppose that a state of a system of A bosons is 



spanned on M orthonormal modes u»(x), so that the 
bosonic field operator can be written as 
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where &j annihilates a boson from the mode Ui{x). The 
joint probability density for A" bosons in the state is 
given by 



(JV) 
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&(x 1 )...W(x N )$(x N )...$(x 1 )\il>). (1) 



This is equivalent to the modulus squared of the A^-body 
wave function. The task is to repeatedly generate a set 
of N numbers £i, according to (JTJ. As noticed by 

Javanainen and Yoo [2|, the probability density above 
can be decomposed into a product of one-dimensional 
conditional probabilities 
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If one desires to generate a set of A positions according to 
the probability pf^P , it is enough to first pick a position 

£i with the probability p^(x\) then a position 



x 2 
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£2 with the pLu^alCi) and so on. 
Each one-dimensional probability density of @ takes 
the following functional form 
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where the coefficients ^ = are calculated from 
with xi = £1, ...,x r -i = £ r _i. 

The important observation 0,0 is that for large N the 

(r) 

coefficients ch assume approximately constant values af- 
ter certain critical number N cr u of positions have been 
drawn, i.e. Cy = consty (r) for r > N cr u. Assuming 
that N cr it <C N, these constant values vary from one po- 
sition measurement of N particles to another. In other 
words, first N cr a particles determine the shape of the 

(r) 

one-dimensional probability densities pYA . The explana- 
tion of this "localization" of values of the coefficients is 
provided in the following. 
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THE CRITICAL NUMBER OF ATOMS. 



particles saturates this inequality 



A condensate is the quantum many-body state of N 
bosons in the form (c') N \0), with denoting a creation 
operator of a particle in a single-particle wave function 
c{x). 

Let us consider a superposition of two condensates 



\ij) = (c\) N \Q) + (cir\0), 

where the operators c\ and £2 annihilate a boson in nor- 
malized single-particle wave functions c\(x) and 02(2), 
respectively. These wave functions need not be orthogo- 
nal. An unimportant normalization factor of \ip) is omit- 
ted. Let us also define a spatial overlap of c\ and C2 



|ci(x)||c2(x)|di. 



(4) 



This quantity assumes values from the interval [0, 1]; it 
vanishes when the wave functions are spatially separated 
and is equal to 1 when |ci(a;)| = |c2(x)|. In the later case, 
a density measurement does not distinguish the conden- 
sates and the corresponding coefficients (j^J would never 
converge to constant values. Bellow we focus on the cases 
where |ci(x)| ^ |ca(ar)|. 

The A-dimensional probability density for this state is 
given by 



JY 
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P\S « ^Mxi)] 2 + l[\c 2 (x i )\ 2 +2]l\c 1 (x i )\\c 2 (x i )\ 
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cos ( ^(arg(ci(a;j)) - arg^fo))) 



(5) 



One sees that the total probability due to the last term in 
(J3J) is limited from above by 2o N , while the first and the 
second term give 2. Thus, if |ci(ac)| 7^ |c 2 (a:) | then o < 1 
and o N — ► when N increases, so the last term can be 
neglected. The remaining two terms define two sectors 
of the A-dimcnsional space where the probability density 
is substantial. Again, if |ci(x)| ^ |c2 (ac) | then these sec- 
tors are spatially separated for large enough N. It can 
be easily seen, by extending the definition of the over- 
lap J3J to multidimensional products of single-particle 
wave functions, JTI^Li l c i( a; i)ll c 2( a; i)|da;i. Such an over- 
lap is given by o N and vanishes exponentially fast with 
N. Finally, if one randomly chooses a point from the N- 
dimensional space according to the probability (JSJ), the 
point will belong to the sector described by either the 
first term or to the sector described by the second term 
in JSJ . Thus, one samples either the condensate in c\ (x) 
or the condensate in C2(x). 

The natural requirement for this separation of the two 
sectors is the small value of the overlap in A-dimensions. 
Namely, o N < 0.01, for example. The critical number of 



In (0.01) 
ln(o) : 
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where the number 0.01 is arbitrary (must be less than 1) 
and due to logarithm in the numerator of JSJ), the N crit 
is not sensitive to its precise value. The (0 provides an 
estimate of an average N cr i t over many repeated mea- 
surements. The critical number of particles for a single 
measurement can differ significantly from (|rj|). 

To better understand the meaning of N cr i t consider 
the case of the Gaussian wave functions 



| Cl (x)| 5 



|c 2 (*)| 5 



The overlap is a function of the parameters L and a 



(7) 



From © and 



N crit = - ln(0.01) 



L 2 ' 



Reversing the problem one can ask, which Gaussians can 
be distinguished given the total number of particles N 
and the width a. The answer: these displaced by more 
than 2a-J — ln(0.01) /N. So, if the state \ip) was an equal 
superposition of all of Gaussians, labeled by the parame- 
ter L, then a single measurement of positions of N parti- 
cles will converge to a subset of Gaussians centered within 
the distance cry — ln(0.01) /N from some random point 
L = L . 

Now, it is clear why the coefficients in (0) assume con- 
stant values only approximately. The "localization" of 
their values improves with the number of already drawn 
particles r like Xj^fr. 

Another important example involves two condensates 
that show shifted interference-like patterns 

IciMI 2 = 5- (l+cos(z+0)), |c 2 (a;)| 2 = -^(l+cos^+p)) 

Z7T Z7T 

for x G [0.271"]. Assuming — 7r < p < (j> < it, the overlap 
of ci and ci is a function of the relative phase 4> — p 



0= 1 
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For small phase shifts <j> — p <C 1 the distinction of the 
two condensates requires at least the number of particles 



81n(0.01) 

{<t>-p? ■ 



This example is important because N cr u above is a crit- 
ical number of atoms for the measurement to "decide" 
which interference pattern is realized in a single density 
measurement of two colliding condensates 0, Q ■ 
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THE METHOD 

The previous part of this work establishes that for 
N ^S> N cr i t almost all the positions of particles in a sin- 
gle simulation are randomly drawn from the same one- 
dimensional probability density. The single measurement 
on the state \ip) can be formally replaced by the measure- 
ment on a Bose-Einstein condensate - the state where 
all the particles are described by the same single-particle 
wave function. Indeed, the measurement of particle po- 
sitions performed on a condensate consists on repeated 
sampling of a single probability density. 
This idea leads to the following postulate: 
The act of position measurement on a many-body state 
\ip) is equivalent to the measurement performed on a con- 
densate with the single-particle wave function 

M 

c{x;Qi, ...,Qm) = £ QiUi(x), (8) 
j=i 

where Qi = qie lBi are random complex parameters satis- 
fying the normalization condition 

M 

£lQi| 2 = l- (9) 

i=l 

The simulation of a single measurement of posi- 
tions of N particles is now reduced either to ran- 
dom generation of N positions, each with the prob- 
ability density \c(x; Q\, Qm)\ 2 , or to usage of the 
N\c(x; Qi, Qm)\ 2 as a smoothened density of mea- 
sured particles. 

The only thing left is to provide a probability distribu- 
tion of random parameters Qi. This distribution arises 
from the projection of the iV-body state of interest on 
a condensate described by the wave function (JHJ). The 
operator that annihilates a particle from the condensate 
mode c(x; Q%, Qm) reads 

M 

c = j2q:^- (io) 

i=l 

The condensate of N particles is given then by 

1^,...^)= * (cty>). (ii) 

Finally, the probability distribution of Qi 's for the state 

p(Qi,... ! Qm)=^>I(c;Qi,-,Qm|V)I 2 , (12) 

with the constrain The normalization factor A/j^,) 
has to be computed for each state \ijj) separately, so that 
/ P{Q x ,...,Qm)AQ x ..AQm = 1. 

The proposed method consist of two stages: (i) random 
generation of M complex parameters Qi from the prob- 
ability distribution (|12f> and under the constrain @, (ii) 



using this set of QiS the wave function (JSJ is constructed 
and its modulus squared is taken either as a probabil- 
ity density for generating N particle positions or as a 
smoothened density of particles in a single measurement. 

At this point three issues need to be clarified. First, the 
condensates <|11[) form an overcomplete and nonorthogo- 
nal basis in the Hilbert space of all iV-body bosonic states 
spanned on the modes Ui(x). The nonorthogonality and 
overcompletness of the condensates are vital if they are 
to describe the outcomes of measurements that are also 
" nonorthogonal" and " overcomplete" , in the sense that 
a weighted sum of two or more outcomes is also a valid 
outcome. Conversely, if the basis formed by the conden- 
sates was orthogonal then every result of measurement 
would be described by one of the basis states and the 
outcomes corresponding to combinations of condensates 
would have been unjustly excluded. 

Second, the postulate above relies on the observation 
that for large number of particles N the system in a state 
\tp) is driven to a condensate state by the position mea- 
surement. In an actual experiment, the positions of all 
particles are measured at the same instant and the ini- 
tial state is destroyed due to absorption of light by parti- 
cles. It is conceivable, though, that starting from a state 
with macroscopically occupied modes one can create a 
condensate by measuring positions of a fraction N cr it of 
particles. In the case of interacting particles, such a con- 
densate might be short lived since it is not an eigenstate 
of the Hamiltonian of the system. 

Third, the form of a possible outcome of a single den- 
sity measurement ||SJ| resembles the wave function of a hy- 
pothetic condensate in the spontaneous symmetry break- 
ing approach, see 0] pages 3 and 4. This work provides 
yet another justification for this approach: It gives cor- 
rect description of measurements if N 3> N cr a. 

EXAMPLES 

Fock state. Lets check how the new method works 
for the Fock state 

|V>) = |ni, ...,n M ). 

From 

P(Q 1 ,...,Q M ) oc IQl 1 x ... x Q n A f\ 2 , (13) 

where M is a number of modes of the state. The proba- 
bility above does not depend on phases 9i of Qi's, so the 
phases will be uniformly and independently drawn from 
the interval [0, 2tt). Having said that, one can restrict 
the probability density to the M -dimensional space of 
moduli qi. The maximum of the probability (under the 
constrain (JHJ) is located at qi = \J ni/N. For large num- 
ber particles per mode the probability is sharply peaked 
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around the maximum point. In fact, one can fix the val- 
ues of moduli qi = \J rii/N and simply draw the phases. 
Thus the expected density of particles found in a single 
experiment is given by 



\c(x;( 




Uj (x) 



This example can be immediately extended to the mix- 
tures of Fock states. 

In particular, when one considers the state \tp) = 
\N/2,N/2) as in Q, with m{x) = ^=e" and u 2 (x) = 



—7==e lx one gets the prediction for a single measurement 



\c( X ; e u e 2 )\ 2 = -L(i + cos(x + {e l - e 2 ))). 

Z7T 



Since both phases are random, their difference is also ran- 
dom and the results from Q are recovered: Each density 
measurement will reveal interference fringes shifted by a 
random phase. 

Schrodinger cat state. Since the presented method 
has ambition to generate outcomes of a single measure- 
ment performed on any state with macroscopically occu- 
pied modes, it should also work for the Schrodinger cat 
state 

\il>) = \N,Q) + %N). 

If corresponding modes u\(x) and u 2 {x) are spatially sep- 
arated, the exact calculation shows (see the section where 
N cr it is derived) , that the single measurement will reveal 
either all particles in the first mode (in the state |iV", 0)) 
or all particles in the second mode (state |0, N)). 
The probability density for Q^s 



P(Qi,Q2)cxqi 



2N 



2q?q? cos(JV(0i 



0), 



under the constrain ©, possesses two well isolated global 
maxima at (qi = l,q 2 = 0) and (qi = 0, q 2 = 1) inde- 
pendent of the phases Qi. For large number of particles 
the probability is appreciable only in the narrow vicin- 
ity of these global maxima. This means that one can 
see either "alive" or "dead" cat in a single measurement 
and that superposition of " dead" and " alive" is strongly 
suppressed. The approximate character of the method 
manifests itself in the fact that even if the modes U\{x) 
and u 2 (x) are spatially separated, there is nonzero prob- 
ability for the superposition. 

Bogoliubov vacuum. A more complex but very im- 
portant example is the Bogoliubov vacuum state. For 
this state, there is also another method of predicting out- 
comes of a single density measurement available . The 
Bogoliubov vacuum is especially useful to describe evolv- 
ing Bose-Einstein condensates. It can be expressed in the 
particle representation as 



\ia\a\ 



t i \ N/2 
+ ^%% |0)> 



where Aj's are real numbers, see [5|,|6||. The corresponding 
probability function assumes the following form 



P(Qi, ...,Qm) oc 



N/2 



vi=l 



In the case of just one depletion mode, M — 2, Ai 
1 , | A2 1 < 1 and the probability becomes 



P({Q}) cx (q$ + \\q\ + 2\ 2 q\q\ 008(2(0! - 2 ))) 



N/2 



The last expression shows a single peak at (q± = 1, q 2 
0) of the variance (Aq 2 ) 2 cx 1/N. 



SUMMARY 

The method presented above consists on sampling of a 
2A/-dimensional probability density (|12|l . What one has 
to draw is a set of M complex numbers Qi to determine 
an outcome of a single measurement. This sampling is 
much simpler, however, than the direct sampling of parti- 
cle positions from an iV-dimensional probability density. 
First of all, for the states with macroscopic occupation of 
modes M <§; N and the dimensionality of the sampling 
problem is strongly reduced. Second, as shown in the 
examples above, the probability distribution l|12(l often 
assumes a simple shape and either analytical calculations 
can be performed or the importance sampling numerical 
techniques can be applied. 

The justification of the method provides also a link 
between sampling of a iV-body probability density and 
the corresponding spontaneous symmetry breaking guess 
for the outcome of the single density measurement. 
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